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ABSTRACT 



It is known that near-Earth objects (NEOs) during their orbital evolution may 
often undergo close approaches to the Sun. Indeed it is estimated that up to ~70% 
of them end their orbital evolution colliding with the Sun. Starting from the present 
orbital properties, it is possible to compute the most likely past evolution for every 
NEO, and to trace its distance from the Sun. We find that a large fraction of the 
population may have experienced in the past frequent close approaches, and thus, as a 
consequence, a considerable Sun-driven heating, not trivially correlated to the present 
orbits. The detailed dynamical behavior, the rotational and the thermal properties of 
NEOs determine the exact amount of the resulting heating due to the Sun. 
In the present paper we discuss the general features of the process, providing estimates 
of the surface temperature reached by NEOs during their evolution. Moreover, we in- 
vestigate the effects of this process on meteor-size bodies, analyzing possible differences 
with the NEO population. We also discuss some possible effects of the heating which 
can be observed through remote sensing by ground-based surveys or space missions. 

Key words: minor planets, asteroids - meteors, meteoroids - Solar system: general 



INTRODUCTION In this paper we discuss a new heating process, which has 

• i- J ■ been neglected so far, namely the heating due to close ap- 

, It is well known that the physical properties of most aster- , . jn „ „, . „ . , , , „ , ,. 

y^Nj , \ rr proaches to the bun. this effect may be relevant tor bodies 

s_J oids (and their fragments, the meteorites) have been affected , , ., ,. ,. , 

£2 such as NEOs, whose penhefia are sometimes very cfose to 
by heating processes. The analysis of the available meteorite 

collections shows that meteorites and their asteroidal par- 
ent bodies have been thermally altered to some degree. The 
alteration ranges from low temperature aqueous processes 



the Sun, and significant temperatures can be attained. 

In the present NEO population the fraction of bodies with 

relatively small perihelion (q) is very small. However, this is 

by no means representative of the past (and future) dynam- 
(<300 K) up to partial or complete melting, differentiation . , , , „ , , . , , . , \ . , „ 

, , — ni ~x icaf tracks. Detailed simulations show that a large traction 

and fractionaf crystallization (>1150 K) (Kcil 2000). The r ,. TT ^ , , , ,, ., , ,. 



of NEOs may have had small perihelion distances for some 

, , . , , , , , , time, hence experiencing episodes of strong heating. More- 

mg induced by decay of short-lived radionuclides, notably . . . , , , ., , , ,. ... 

26 a HT~- o — d II u . „ . over, since meteoroids have an orbital evolution similar to 



main source of heating has been identified as internal heat- 



3 A1 |Grimm fc McSween|[l993l ). The effects of heating pro- 
cesses have also been identified on the surfaces of asteroids 
by remote sensing spectroscopy of primitive main belt aster- 



that of NEOs, this process may also affect the meteorite 
samples. 

, . , .. , The detailed effects of the heating during Sun close ap- 

oids (tiiroi ct al. 199o; Vilas and Gaffcy 1989). Despite the , , , , „ , ,. 

j- T~j Z — ■ i proaches depend on the thermal properties of the bodies. 



general consensus on radionuclides heating, other processes 
have been suggested to occur: for instance heating during 



Albedo, emissivity, macroscopic surface roughness and ther- 
mal inertia are the physical parameters that determine the 
the preaccretion phase or due to collisions. The latter may r „ , , 

f , . , . , surface temperature. 1 he fatter parameter measures the re- 

produce high temperatures m the regions involved m high . . „ . , , . , , 

sistance of a materiaf to temperature changes, e.g. due to 
energy impacts, often causing the partial melting of the in- . , / . , . .„ . ,. , , ... 

| ^ r , — — -i| n i the varying day/night illumination caused by the rotation 

volved bodies (see discussion m Rivkm et al. 20021). „ , , ™ ... r , , , „ , ., 

v ' of the body. Ihe rotation state of the body afso affects its 

temperature profile. For instance, for given set of thermo- 
physical parameters, a rapid rotating asteroid has a tem- 
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perature distribution more smoothed out in longitude than 
a slower rotating one. The direction of the body's spin vec- 
tor causes also seasonal effects. 

In this work we estimate the amount of heating due to the 
Sun that NEOs and meteorites are likely to have experi- 
enced. We do this by studying their dynamical history. We 
then use thermal models to correlate their dynamical his- 
tory with surface temperatures. Finally, some possible con- 
sequences of the heating are discussed. 



2 ORBITAL PROPERTIES OF NEOS 

NEOs evolve rapidly in the Keplerian orbital elements space, 
(a,e,i), due to a combination of close encounters with the 
terrestrial planets and resonances with the giant planets. 
Two aspects of this orbital evolution are relevant for the 
analysis of the heating process: (i) The probability that each 
NEO attained a q below a given limit (q s ) during the evolu- 
tion; (ii) The cumulative time spent at distances below q s - 
In the following sections, we first explain how we estimated 
these quantities and then present and discuss the results 
that we have obtained. 

2.1 Estimating the past dynamical histories of 
NEOs 

The most complete mode l of the orbi t al distribut ion of 
NEOs is that developed in iBottke et~aT1 (|2000l . I2002T ). This 
model assumes that NEOs come primarily from 5 intermedi- 
ate sources: the 3:1 mean-motion resonance with Jupiter, the 
ve secular resonance, the Mars-crossing population (IMC), 
the collection of resonances crossing the outer asteroid belt 
(OB) and the trans-Neptunian region (as dormant Jupiter 
family comets; JFC). In the NEO orbital space (q < 1.3 AU), 
the steady state orbital distribution of the bodies coming 
from each of these sources (P S ourcc(a, e, i)) was computed 
using a large number of numerical simulations. The orbital 
distribution of the NEO population was then constructed 
as a linear combination of the distributions related to each 
source, that is: 

i?NEo(a, e, i) = N s 

ource Psource (fl, 6, i) . 

source 

The coefficients iVsourcc for this combination were deter- 
mined by fitting the distribution of the NEOs discovered or 
accidentally recovered by the Spacewatch survey, once ob- 
servational biases had been taken into account. This model 
was later shown to reproduce adequately also the orbital dis- 
tributions of NEOs discovered by the LINEAR and Catalina 
Sky surveys (|Zavodnv et alj|2008h . 

Given a specific NEO, using the Bottke et al. model we 
estimate the probability that it comes from each of the con- 
sidered intermediate sources, as follows. We first define a box 
in (a, e, i) orbital space in which the NEO currently resides. 
Then, for each source, we compute Psourcc, as the integral 
of -Rsourcc {a, e, i) over the box; we also compute Pneo in 
an equivalent manner. Then, the probability that the NEO 
comes from a given source is given by: 

P(source) = N BOUICC R BOUTCC / R^eo ■ 
In our calculation, we use boxes of size 0.1AUx0.1x5°. The 



probability P{q s ) that the considered NEO achieved in the 
past an orbit q < q s is 

P(q s ) = /J P(source)P S ourcc(g s ) , 

source 

where P S ource(g) is the probability that objects coming from 
the source achieved q < q s before entering the (a,e,i) box. 
Similarly, the mean time T(q s ) spent on orbits with q < q s 
can be computed as 

T(q s ) = P(source) 

1 source 

source 

where T SOUICC (q a ) is the mean time spent on orbits with 
q < q 3 by the particles reaching this kind of orbits from 
the source, before entering the (a, e, i) box. To compute 
Psource (q s ) and T S ource(<Zs) we went back to the original simu- 
lations used by Bottke et al. to construct the Psource (a, e, i) 
distributions. The computation is just a matter of book- 
keeping, while the outputs of these simulations are read se- 
quentially in time. For simplicity (and availability of the 
original simulations) we did this only for the 3:1, va and OB 
sources. We then assumed that the objects coming from the 
IMC population behave statistically as a 1:1 combination of 
the objects coming from the 3:1 and vq sources; similarly 
we assumed that the objects coming from the JFC source 
share the same statistical properties of those coming from 
the OB source. These assumptions are of course not exact, 
but are relatively close to reality. By doing this, the equa- 
tions defining P(q s ) and T(q s ) are effectively restricted to 
three sources (3:1, v&, OB) with new coefficients 

P'(3:l) = P(3 : 1) + 1/2P(IMC) , 

P> 6 ) = P(j/ 6 ) + 1/2P(IMC) , 
P'(OB) = P(OB) +P(JFC) . (1) 

2.2 Results 

In Figure [T] we report the a — e scatter plot of the NEO 
population for several values of q„ . For each panel, the color 
of NEOs is coded according to the probability of experi- 
encing q < q s . The curve q — q s is also overplotted. The 
figure clearly shows that the orbital paths followed by each 
NEO very often have led to small perihelion distances, much 
shorter than their present q. In particular, this is evident 
for q s = 0.05, 0.1 AU: only few NEOs are presently below 
those values, but a significant fraction of the NEO popula- 
tion spent some time, during the evolution, with a smaller 
perihelion. The average probability vs q s is shown in fig. [5] 
(left panel). 

It is obvious that the present orbit and the past history are 
correlated (this is the rationale of our computations); how- 
ever, the correlation is not trivial. As a general trend, P{q s ) 
decreases for objects with increasing value of their current 
perihelion distance q. Only bodies with semi-major axes in 
the region of the main resonances v§ and 3:1 (a = 2—2.5 AU) 
can have large P(q s ) even if they currently have a large value 
of q (see fig. [TJ . This happens because these resonances en- 
tail large oscillations of the eccentricity. Moreover, we re- 
mark that the bodies with the smallest value of P(q 3 ) are 
those with large semi-major axis (beyond 2.5 AU), mostly 
coming from the OB source. We may thus expect that a sig- 
nificant percentage of NEOs, originating in the inner part of 
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the Main Belt and delivered through the v§ resonance, have 
been affected by the consequences of the heating process, 
and thus possibly physically altered (see also the discus- 
sion in the following). The other bodies should have kept 
more systematically their original properties, at least for 
what concerns the alterations due to heating. Another in- 
teresting result concerns the time spent at distances smaller 
than q 3 . In Figure [3] we show the NEO population with 
colors coded according to the cumulative time elapsed for 
q < q s . We found that, even for small q s , the time elapsed 
by NEOs at close distances to the Sun can be considerably 
high, reaching the 10% of the typical life times (10 Myr) 
or, in a few cases, even more. The average cumulative time 
spent below q s is shown in fig. [2] (right panel) . 
A closer look at the results of our simulations in order to 
identify the most heated NEOs shows that 1% of NEOs have 
a 50% probability to have been at q < 0.1 AU. Such proba- 
bility becomes 5% at 0.2 AU. Their cumulative time spent 
below such distances can vary considerably, from several 
10 3 yr to several Myr. Interestingly, one of the object pre- 
dicted to be highly processed is the peculiar 3200 Phaethon, 
whose probability to have been at q < 0.1 AU is 56% for a 
cumulative time of 0.3 Myr. 

On the other hand, there are NEOs that never went close 
to the Sun, and therefore they represent the least processed 
bodies. For instance, 2% of NEOs have a probability of less 
than 10% to have spent time below q s — 0.8 AU. These two 
extreme sample of NEOs -the hot and cool bodies- show 
interesting features in the a — e plot (see fig. 13}. 



3 THERMAL PROPERTIES OF NEOS 

In the previous section we have shown that a fraction of the 
NEO population has experienced small perihelion distances; 
some of these bodies had q < q s also for a long cumulative 
time. For these objects the effects of the solar heating may 
have altered their surface properties. In this section we use 
asteroid thermal and thermophysical models to investigate 
in more detail the dependence of the surface temperature 
of a body on its heliocentric distance and thermophysical 
properties. Obviously the presence of surface regolith, and 
its continuous mixing , due to microcollisions or tidal mixing 
dMarchi et al.ll2006al ) or the YORP reshaping (|Harris et alJ 
l2009h . may -at least partially- mask the above mentioned 
alterations. 

The surface temperature depends on the body's thermal and 
rotational properties. However, rough estimate of surface 
temperatures at perihelion can be obtained from the equa- 
tion of instantaneous thermal equilibrium with sunlight: 

r=[(i-^ 0g -VVV- 1 ] - 25 (2) 

where A is the bolometric Bond albedo, Sq is the solar con- 
stant at 1 AU (1329 Wm~ 2 ), e is the infrared emissiv- 
ity, a the Stefan-Boltz mann constant and n the so-called 
beaming parameter (see lHarris and Lagerros" |2002| . and ref- 
erences therein). The latter can be seen as a measure of 
the departure of the asteroid temperature distribution from 
that of a spherical, smooth body with all surface points 
in instantaneous thermal equilibrium with sunlight (which 
would have r\ — 1). The value of r\ is a strong function of the 
macroscopic roughness and the thermal inertia of the surface 



i|Spencer et al. II 19891 ; lHarris lll998l ; ll5elbo' et al. l2007h . The 
latter parameter, which measures the resistance of a mate- 
rial to temperature changes, is defined as F = ^/ptc, where 
r is the thermal inertia, p is the density of the material, c 
is its specific heat content and k the thermal conductivity. 
In particular, r\ > 1 indicates a surface with a temperature 
cooler than the one of the instantaneous thermal equilib- 
rium, whereas rj < 1 is expected for a surface with low ther- 
mal inertia and significant roughness. 

An upper limit for the NEO surface temperature can be ob- 
tained by setting in Eq. ([2]) a very low bolometric albedo 
(A=0.01, corresponding to a geometric visibl e albedo, py, 
of 0.03 and the default phase integral -G=0.15. lBowell et al.l 
(1989)-; pv = 0.03 is at present the lowest geometric albedo 
measured for NEAs), a ze ro thermal inertia a nd a strong 
infrared beaming f 17=0.57. IWolters et al. 1120051 ). For an in- 
frared emissivity of 0.9 -commonly adopted for the min- 
eralogy of NEOs ilLim et al. 1 120051 ; ISalisburv et al. I Il99ll ; 



iMustard and Hays I Il997l ) 



we find temperatures of 1968 
and 1391 K at 0.05 AU and 0.10 AU respectively. Aver- 
age parameter values for NEOs (i.e. A=0.05 from pv=0.14 
IStuart and Binzel 112003 ) and G=0.15; t) = 1.0) yields sub- 
solar temperatures of 1782 and 1260 K at 0.05 AU and 
0.10 AU, respectively. A lower limit for the surface temper- 
ature can be obtained for a small body with instantaneous 
re-distribution of the heat in the whole volume, with no dif- 
ference between the temperature of day and night: namely 
an isothermal object. This latter case can be obtained using 
77=4 in Eq. J2j. We obtain 1343 K at 0.05 AU and 949 K at 
0.10 AU for yl=0.05. We note that the la rgest ever measured 
value for a NEO is 3.1 (3671 Dionysus; H arris and Davies I 
1999). Figure \5\ shows the surface temperature as function 
of the heliocentric distances for the three cases described 
above. 

However, observable consequences of heating can come out 
if and only if the whole surface (or a large amount of it) has 
been involved, and thus altered. In order to estimate the 
fraction of the surface above a certain temperature thresh- 
old, more accurate knowledge of the temperature distribu- 
tion on the surface of asteroids is needed. Thermophysical 
modeling is required to obtai n this informatio n. A ther- 



Spencer et al. 1 1 198 



mophysical model (TPM. see lH arris and Lagerros 



Lagerros 



19961 ; 



Emery et al. 





2002 




1998 



D- ' elbo' and Tan ' ga I T^rT^ef^ce s therein) de scTibts 



an asteroid as a polyhedron made by a mesh of planar facets. 
The temperature of each facet is determined by numerically 
solving the one-dimensional heat diffusion equation into the 
subsurface with boundary conditions given by the diurnal 
variable illumination and by the energy irradiated away by 
each facet at the surface, and by setting the heat flow at 
the deepest subsurface element equal to zero. In our imple- 
mentation of the TPM, the subsurface is divided into 32 
slabs of thic kness 0.25 tim es the diurnal heat penetration 
depth (l s = \/ fi;p _1 c _1 aj _1 ), where u) = 2ir/P with P the ro- 
tation period of the asteroid; the deepest subsurface element 
is therefore 8 l s beneath the surface. Surface and subsurface 
temperatures are controlled by a number of physical param- 
eters including the heliocentric distance q, A, e, T, P, the 
direction of the rotation axis with respect to the Sun, and 
the shape of the body. 

In the following, we assume a spherical asteroid with the 
spin vector perpendicular to its orbital plane, P=6 hours, 
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e=0.9, A=Q.l, T=200 Jm~ 2 s~ 05 K~ 1 . The latter is the 
average value of th e thermal inertia for km-sized NEAs 
ijDelbo' et al. II2007T ). The TPM was run on circular orbits 
with heliocentric distances ranging from 0.05 AU to 0.5 AU. 
At the beginning of each run, a complete rotation of the 
body is performed and the average temperature of each 
facet is recorded. The temperature at each rotation step 
(usually 360 steps per rotation are performed) is calculated 
assuming instantaneous thermal equilibrium between ther- 
mal infrared emission and absorption of solar energy, (i.e. 
Ti = [(1 - A)S e q~ 2 cosOiitfe-'ta- 1 ] ' 25 ), where Ti is the 
temperature of the i t h surface element of the mesh, q is the 
heliocentric distance of the body, and 9i is the angle formed 
by the normal of the i t h surface element of the mesh with 
the direction to the Sun, and it varies with the time t. The 
average temperature of each facet is then used as the initial 
temperature of all 32 subsurface slabs including the one at 
the surface. The temperature profile as a function of time, 
i.e. the asteroid rotational phase, is monitored on each facet 
during the warm up phase, which can take up to 50-100 (de- 
pending on the value of q) full rotations until the tempera- 
ture profile stabilizes and the initial temperature conditions 
are forgotten. A final rotation is then performed and the 
maximum temperature of each facet is stored. The areas of 
those facets whose maximum temperatures are found above 
a given temperature threshold (T„) are added up and the 
ratio of this area to the total area of the body surface is 
calculated for a range of values of T s from up to 2000K. 
The value of T s for which the ratio of the surface whose 
maximum temperature is above T s is equal to 0.7 and 0.5 
is found by interpolation. Figure [5] shows these values of 
T s for different heliocentric distances: these values represent 
the temperature above which a 70% and 50% fraction of the 
surface area of the body was heated to. 
We also note that the temperature above were calculated us- 
ing a TPM with smooth surface. However, it is well known 
that roughness increases the average surface temperatures 
by a maximum of 20-30 % with respect to those of a smooth 
surface. Running the TPM on orbits with substantial eccen- 
tricity does not affect significantly the final result. 



4 DIFFERENCE BETWEEN NEOS AND 
METEORITES 

Meteorites are collected on the Earth surface after having 
passed some time in a NEO-like orbit. Their history may 
be various: some may be relatively young fragments from 
cratering or catastrophic collisions, or from tidal shattering 
of NEOs, others may have a longer independent history, 
maybe a long orbital evolution as NEOs. Despite of the 
different size, observed NEOs and meteorites are believed 
to share a common origin; moreover the orbital evolution of 
meteorites is typically dominated by the same dynamical 
processes as NEOs. 

However, for what concerns the alterations due to heating, 
we expect some important differences with respect to 
NEOs: (i) Meteorites might have spent in near-Earth space 
a shorter time than the average NEO, because of their short 
collisional lifetime. This is suggested by the fact that the 
statistics of falls times (in the afternoon vs. in the morning) 
is skewed towards afternoon falls with respect to that of 



the expected NEO impacts l|Morbidelli fc Gladmanlll998f ). 
As a result the time spent might have been not enough for 
considerable heating; (ii) Meteorites are systematically by 
far smaller than NEOs. Thus also the thermal history is 
likely different; (iii) bodies that collides with the Earth (like 
meteorites) have preferentially q ~ 1 AU. This introduces a 
further bias with respect to the observed NEO population, 
where many of the bodies have a very low probability of 
collision with the Earth. Indeed, since P{q B ) decreases with 
increasing q, this implies that meteorites are likely less 
heated than NEOs. 

In the limit of a very small body we can imagine an almost 
instantaneous re-distribution of the heat in the whole 
volume, no day-night difference, a quasi-isothermal object. 
If we use the same formulas and basic assumptions to 
explore this extreme case, the factor 4 between the surface 
of a spherical body and its cross section introduces a factor 
4 -1,/4 in the attainable temperature, which consequently 
decreases by ~30% (see Eq. [2] and the discussion in the 
previous Section). Real meteorites may behave in a more 
complex way, but somehow similar to the above scenario. 
In fact, an accurate thermal model suggests that a thin 
layer, in the region close to the subsolar point, may attain 
large temperatures, similar to those estimated for asteroids. 
However, we expect that those thin layers are systematically 
erased during the approach of the meteorite to the Earth, 
mainly due to atmospheric friction. Thus the meteorite "as 
we see it", i.e. after we recover it on the ground, retains 
no effect of the strong surface heating, and its thermal 
evolution can be analyzed in terms of the "isothermal" 
approximation. Thus the thermal history of NEO and 
meteorites can be considerably different. We must take this 
consideration into account when attempting to compare 
laboratory measured quantities on meteorites with analogue 
quantities derived from remote sensing of asteroids. 



5 DISCUSSION: CONSEQUENCES OF THE 
HEATING 

The purpose of the present paper is to point out the 
importance of the heating history of NEOs and meteoroids, 
due to very close passes to the Sun. Even if the general 
ideas seem rather robust, the physical consequences have 
to be explored in deeper detail. In fact a number of 
interesting phenomena can be caused by the mentioned 
heating. Generally speaking the heating can drive volatile 
release, therefore the surface composition of NEOs and, to 
a lesser extend, meteoroids become progressively depleted 
of volatiles. In the case of temperature higher than 1200 K, 
silicate sublimation is expected. 

For temperatures below the melting points of the different 
silicates known to be present on asteroid surfaces, it is worth 
to note that thermal heating can cause annealing: this is 
in general a transition from an amorphous to a crystalline 
structure. The total time that an object spent with q < q s 
is extremely important in this kind of process. This is 
because the rate, k, at which a material is transformed 
from amorphous to crystalline is function of temperature 
via an equation of the form: k = aexp(—E a /bT), where o 
and b are constants determined from experiments (see e.g. 
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iBrucato et al.| [2003. and references therein) and E a is the 
activation energy for the transition. The latter depends on 
the chemical nature of the silicate. 

Concerning the most "primitive" NEOs, i.e. those be- 
longing to the C-types, more severe effects are foreseen. 
For instance the aqueous altered carbonaceous chondrites 
are strongly affe cted by moderately high temperatures 
ijHiroi et al.lll996T ). We note that the percentage of C-type 
NEOs with aqueous alteration features (in particular the 
0.7 fim feature) is rather high among MB As, of about 50% 
|Bus fc Binzelll2002l ; ICarvano et alj|2003h . This percentage 
drops below 10% for NEOs: only 3 objects over more then 
30 C-type NEOs having visible spectroscopic measurements 
show subtle absorption fea tures a t 0.7 /xm. These results 
have been also found by IVilas I (|2005h on the basis of 
ECAS photometry. On the other hand, the hydration is a 
common feature in the meteorite collection, in agreement 
to what our considerations predict. However, the data are 
too few and probably not unambiguous enough to support 
the explanation of the difference of the NEOs hydration 
properties solely in terms of heating. For instance, one 
(2002 NX18) of the C-type NEOs with hydration features is 
probably coming from the Outer Belt, thus has presumably 
been poorly exposed to heating. A similar result holds also 
for 162173, while the dynamical history of 2002 DH2 is less 
clear. An alternative explanation can be that NEOs sample 
peculiar source regions. 

Similar considerations apply also to the cometary objects, 
namely those NEOs coming from the JFC source region 
(about 70 NEOs are estimated to originate in this way). 
The very close passages to the Sun and high temperature 
heating should have produced a rapid release of all volatiles, 
leading to the stage of depleted cometary nuclei. 
Another interesting aspect is related to the space weath- 
ering, which is expec ted to be much more efficient in 
proxi mity to the Sun l|Marchi et al.l l2006bl ; iPaolicchi et al.l 
12007ft . Changing from 1 AU to 0.1 AU increases the 
space weathering rate by two order of magnitude. Similar 
considerations might also apply to other effects depending 
on the solar radiation, such as Yarkovsky and YORP. 
Moreover, we find that the thermal histories of NEOs 
and meteorites may be considerably different, therefore 
the above mentioned processes may introduce systematic 
differences in their physical properties. A future analysis 
of these points will be worthwhile; in general we remark 
that the effects of close passages to the Sun, in terms of 
heating, but also as concerns other processes, should not be 
neglected, especially when evaluating the targets for future 
space missions to NEOs, such as the European/ Japanese 
Marco Polo mission or the USA Osiris mission. 
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Figure 1. Semimajor axis vs eccentricity scatter plot of the present NEO population for several values of q B - NEOs are plotted with 
colors coded accordingly to their probability of stay below q s (white dots indicate zero probability). For each panel, the corresponding 
q = q B curve is also shown. 
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Figure 2. Left panel: average NEO probability of q < q s . The curves for the three source regions have been normalized to the number 
of NEOs coming from each region. Right panel: average NEO cumulative time spent at q < q s . The curves for the three main source 
regions are also shown. 
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Figure 3. Same as FigfT] with colors coded according to the logarithm of cumulative time spent below q s (white dots indicates T s = 0). 
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Figure 4. The disposition of the "hot" and "cool" NEOs. The hot NEOs (red circles) are defined as those having a probability of at 
least 50% to fall below 0.1 AU. The "cool" (blue triangles) are those having a probability of less than 10% to fall below 0.8 AU. 
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Figure 5. Temperature of NEOs as function of their heliocentric distances. In particular: continuous line: temperature of the sub-solar 
point for a flat surface in instantaneous thermal equilibrium with sun light (physical parameters are A=0.05, r\ = 1.0 and e = 0.9). 
Dashed-dotted line: temperature of the sub-solar point for a rough, low-albedo surface in instantaneous thermal equilibrium with sun 
light (physical parameters are A=0.01, r\ = 0.6 and e = 0.9). Dashed line: temperature of an isothermal smooth-surface body (physical 
parameters are A=0.05 r\ = 4.0, e = 0.9). This curve can represent the tempuerature of the meteorites. Dotted lines with symbols: 
temperature above which a 70% (diamonds) and 50% (triangles) fraction of the surface area of an NEO is heated to. See section [3] for 
further details. 



